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ABSTRACT 

This  thesis  reports  the  results  of  an  experimental  investigation  into 
the  spanwise  correlation  of  the  fluctuating  longitudinal  velocity  component 
in  the  early  wake  of  a  flat  plate  placed  parallel  to  a  uniform  stream  and 
allowed  to  vibrate  about  a  vertical  axis  through  the  leading  edge.  The 
experimental  program  provides  for  the  measurement  of  fluctuating  total  head 
signals  at  one  fixed  and  one  variable  spanwise  position  on  the  wake  center- 
line  one  body  thickness  downstream  of  the  trailing  edge.  Analyzation  of 
these  signals  gives  the  spanwise  correlation  of  the  wake  structure  expressed 
as  a  complex  number  and  the  power  spectral  density  of  the  longitudinal 
velocity  fluctuation  along  the  span.  In  addition  the  mean  local  wake 
velocity  is  measured  for  each  value  of  the  ambient  flow  velocity. 

Ambient  velocities  of  8,  11,  and  14.  fps  and  a  single  plate  of  £"  thick- 
ness and  2"  chord,  were  used  for  all  the  tests  which  were  carried  out  in  the 
7^"  x  9"  test  section  of  an  open-circuit  water  tunnel  in  the  Hydrodynamics 
Laboratory  at  MIT. 

An  equation  is  developed  which  relates  the  vortex- induced  forcing 
moment  on  the  foil  to  the  measured  structure  of  the  early  wake.  Using  this 
model,  forcing  moment  coefficients  are  evaluated  as  a  function  of  a  "self 
excitation"  parameter  consisting  of  the  ratio  of  the  mean  transverse  velocity 
of  the  trailing  edge  to  the  free  stream  velocity.  The  role  of  spanwise 
coherence  of  the  wake  is  examined  with  respect  to  its  effect  on  the  coupling 
between  body  vibration  and  forcing  moment. 

Major  Conclusion  of  this  study  is; 

The  increase  in  the  forcing  moment  coefficient,  m.  ,  accompanying  an 

o 
increase  in  the  self  excitation  parameter  is  due  to  increasing  correlation 

or  coherency  of  the  wake  along  the  span  rather  than  to  increasing  strength 

of  the  individual  vortices  being  shed. 
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LIST  OF  SYMBOLS 

2  b  plate  chord;  2  inches 

f  frequency 

f(t)  function  of  time 

FjCnj)    Fourier  transform  of  periodic  signal  of  discrete  frequency  n1 

F2(n2)    Fourier  transform  of  periodic  signal  of  discrete  frequency  n2 
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H  spanwise  length  of  plate;  in. 

H.  fluctuating  total  head;  ft.  of  water 
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i         v^T 

I.  total  moment  of  inertia;  lbs-in-sec2 

K,  total  spring  coefficient;  lbs-in 

m.  average  forcing  moment  coefficient;  dimensionless 

mf  average  amplitude  of  forcing  moment  coefficient;  dimensionless 
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M  vortex  induced  forcing  moment;  lbs-in 

n  discrete  frequency 

N  lift  force  on  plate;  lbs 

p'  fluctuating  static  pressure;  lbs/ft2 

p  mean  static  pressure;  lbs/ft 

t  plate  thickness;  £  inches 

t  time 

u*  fluctuating  longitudinal  velocity  component;  fps 

u  mean  local  wake  velocity;  fps 

V  mean  ambient  velocity  outside  the  wake;  fps 

V  velocity  of  separation  streamline;  fps 


vf  transverse  fluctuating  velocity;  fps 

v  mean  transverse  velocity  in  the  wake;  fps 

V*  — -^  ,  "self  excitation"  parameter;  dimensionless 

w*  spanwise  fluctuating  velocity  in  the  wake;  fps 

w  mean  spanwise  velocity  in  the  wake;  fps 

x  longitudinal  coordinate;  origin  at  trailing  edge 

y  transverse  coordinate;  origin  at  trailing  edge 

z  spanwise  coordinate;  origin  at  mid-span 
angular  plate  displacement;  radians 

o  RMS  of  vibrational  amplitude;  radians 

£  weight  density  of  fluid;  lbs/ft^ 

vortex  circulation;  ft2/sec 

8  phase  angle  of  correlation 

<J>  phase  angle  between  forcing  moment  and  plate  motion;  radians 

Jaa  cross  power  density;  signal  A,  signal  A 

Q  cross  power  density;  signal  B,  signal  B 

^\  cross  power  density;  signal  A,  signal  B 
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I  INTRODUCTION 

A.  Flow- Induced  Vibrations 

Structures  subjected  to  a  flow  field  may,  under  certain  circumstances, 
experience  flow- induced  vibrations.  These  vibrations  are  undesirable  and 
indeed  may  cause  complete  collapse  of  the  structure  due  to  fatigue  induced 
by  the  vibrations. 

Unstreamlined  or  "bluff"  bodies  in  a  flow  field  shed  discrete  vortices 
in  the  wake.  The  shedding  of  these  vortices  gives  rise  to  periodic  lift 
forces  in  a  direction  transverse  to  the  flow  direction.  It  has  been  shown 
in  reference  (2)  that  for  bodies  of  finite  thickness  this  vortex- induced 
lift  force  increases  linearly  with  the  transverse  vibrational  velocity, 
thus  producing  what  may  be  called  a  "self -excited"  vibration.  The  fluid 
mechanism  involved  in  this  feedback  between  motion  and  forcing  force  is 
unknown  but  may  be  expected  to  be  due  to  an  increase  in  the  strength  of 
the  discrete  vortices  being  shed  and/or  to  an  increase  in  the  spanwise 
correlation  of  the  vortices  being  shed. 

B.  Scope  of  the  Investigation 

The  purpose  of  this  investigation  is  to  determine  the  cause  of  the 
amplitude  dependence  of  the  forcing  moment  so  that  engineers  may  devise 
steps  to  prevent  the  buildup  of  large  amplitude  vibrations. 


II     EXPERIMENTAL  EQUIPMENT 

A.  Water  Tunnel 

The  experiment  reported  herein  was  carried  out  in  the  closed-jet,  open 
circuit  water  tunnel  of  the  Hydrodynamics  Laboratory  at  MIT.  The  water 
tunnel  is  fully  described  in  reference  (1).  Figure  1  shows  the  water  tunnel 
test  section.   Pump  number  III  (rated  capacity  3000  G.P.M.  at  a  head  of  45 
ft.  and  1175  R.P.M.)  was  used. 

Flow  measurement  was  obtained  using  the  tunnel  contraction  as  a  flow 
meter.  Its  calibration  was  checked  and  found  to  be  the  same  as  reported  in 
reference  (2). 

B.  Test  Plate 

The  test  plate  used  was  plate  D  from  reference  (2),  having  the 
geometry  shown  in  Figure  15*  The  plate  was  mounted  vertically  in  the  test 
section  in  such  a  way  as  to  permit  rotation  about  the  center  of  its  semi- 
circular leading  edge.  The  plate  is  fabricated  of  aluminum. 

C.  Plate  Motion  Measurement 

1.  Strain  Gages 

Tatnall  type  C9-121-R2C  strain  gages  were  cemented  to  the  torsion 
spring  which  restrains  the  plate  motion  and  were  waterproofed.  The  four 
gages  form  a  wheatstone  bridge  which  may  be  balanced  by  means  of  an  external 
variable  resistance.  The  bridge  output  voltage  is  proportional  to  the  angle 
of  plate  rotation.  See  Appendix  C  of  reference  (2)  for  theory  of  calibra- 
tion of  the  strain  gage  bridge. 

2.  Amplifier 

A  decade  amplifier  model  220-A  made  by  Ballantine  Laboratories  Inc.  was 
used  to  amplify  the  strain  gage  bridge  output. 


3«  Oscilloscope 

An  Eico  model  4.60  DC-wide  band  oscilloscope  was  used  for  visual 
observation  of  the  strain  gage  bridge  output. 

D.  Total  Head  Fluctuation  Measurements 

1.  Total  Head  Turbulence  Gages 

The  gages  for  measuring  the  total  head  fluctuations  (turbulence)  are 
fully  described  in  references  (3  and  4).  Figures  3  and  4,  Appendix  A,  are 
photographs  of  the  total  head  gages  and  Figure  15,  Appendix  A,  is  a  drawing 
of  the  probes  and  test  plate. 

2.  Electrometer  Preamplifier 

Two  electrometers  (one  for  each  turbulence  probe)  were  used  to  match  the 
probe  and  amplifier  impedences. 

The  electrometers  are  fully  described  in  references  (3  and  4)» 
Figure  16  is  a  block  diagram  showing  all  electronic  components. 

3.  Amplifiers 

Two  Dana  model  2000  solid  state,  wideband,  differential,  low-level  D.C. 
amplifiers  with  integral  power  supplies  were  used  to  amplify  the  turbulence 
probe  signals.  The  amplifier  has  a  flat  response  to  -3db  from  DC  to 
approximately  70  KC. 

4.  Oscilloscope 

A  Dumont  type  333  dual  beam  oscilloscope  was  used  for  visual  observa- 
tions of  the  two  turbulence  probe  signals. 

5.  Portable  Tape  Recorder 

A  Precision  Instrument  Co.  portable  tape  recorder  model  PS  200  A  was 
used  to  record  the  turbulence  probe  signals. 


6.  Endless  Loop  Tape  Recorder 

A  Sanborn  model  2000  endless  loop  tape  recorder  was  used  to  give  a 
continuous  input  of  the  turbulence  probe  signals  into  the  wave  analysis 
system. 

7.  Wave  Analysis  System 

To  analyze  the  turbulence  probe  signals,  a  Technical  Products  Co. 
TP  625  Wave  Analyzer  System  was  used.  This  system  consists  of  the  following: 

a.  One  TP  626  Oscillator 

This  oscillator  produces  an  output  frequency  which  is  continuously 
variable  from  97  KC  to  122  KC  over  three  different  frequency  bands: 
97.0  KC  -  97.25  KC;  97.0  KC  -  99.5  KC;  and  97.0  KC  -  122.0  KC.  The 
oscillator  output  is  used  to  modulate  the  signal  to  be  analyzed  in 
the  TP  627  analyzers. 

b.  Two  TP  627  Analyzers 

These  analyzers  determine  the  frequency  and  amplitude  of  the 
individual  components  of  a  complex  wave  within  the  frequency  range 
of  2-25,000  cycles.  The  results  are  measured  in  percent  of  total 
signal. 

These  analyzers  are  equipped  with  two  matched  TP-218  type 
filters  as  required  for  operation  in  conjunction  with  a  TP-645 
Multiplier. 

c.  One  TP-6/*5  Multiplier 

This  multiplier,  when  used  in  conjunction  with  the  other  components 
of  the  system,  gives  the  cross  power  spectrum  of  two  complex  waves  with- 
in the  frequency  range  of  2-25,000  cycles. 


5. 

Ill  THEORY 

A.  Forcing  Moment  Coefficient 

In  reference  (2)  Noutsopoulos  develops  the  equation  of  motion  of  a 
single  degree  of  freedom  plate  spring  system  vibrating  about  an  axis  through 
the  leading  edge  of  the  plate.  The  motion  is  caused  by  a  forcing  moment 
resulting  from  flow  separation  at  the  trailing  edge  of  the  plate.  Assuming 
two-dimensional  motion,  the  equation  of  motion  of  the  plate  may  be  written 

I.CI  +  Ct<*  +  Ktot  =  M 
where  I.  =  the  total  moment  of  inertia 

C.  =  the  total  damping  coefficient 
K.  =  the  total  stiffness  of  the  system 
and  o(    =  the  angular  plate  displacement  in  radians    ex  =  c*^  £ 

From  purely  dimensional  considerations,  the  simplest  form  of  the 
forcing  moment  can  be  written 

M  =  mf  £pH(2b)2V2 

where  m_  is  the  spanwise  average  forcing  moment  coefficient,  H  is  the  span 
and  V  the  ambient  fluid  velocity.  When  the  forcing  moment  is  periodic  as 
in  the  case  of  excitation  by  vortex  shedding,  the  moment  coefficient  is 
conveniently  written  in  complex  form  as 

i(wt  +  4>  ) 

m.  =  mf  e  v     ~  ' 

o 

where  w  =  angular  frequency  of  vortex  shedding  from  one  side  of  the  plate 
in  cps 

4>  =  phase  angle  tetween  the  spanwise  average  forcing  moment  and  the 
motion  of. 

t  =  time  in  seconds. 


6. 

The  moment  on  the  plate  may  also  be  written  in  terms  of  a  lift  force 
and  an  appropriate  moment  arm  where  the  lift  force  is  a  function  of  the 
circulation  which  is  a  measure  of  the  vortex  strength. 

Therefore 

M  ™  2bN 

where  2b  is  the  chord  length  and 
N  is  the  lift  force . 

The  lift  force  may  be  written  as 

N  =  p  VpH 

where  P  is  the  circulation  on  the  plate  and  is  a  function  of  position  along 

the  span  and  time  such  that: 

H 
N  =  p  V  I  T  (z,t)dz  . 
0 

If  the  circulation  is  considered  to  be  a  periodic  function  then 

Pu,t)  =  r0w.ib*  +  * +  9(z)) 

with  U  (z)  being  the  phase  angle  of  the  circulation  at  some  position,  z,  along 
the  span,  measured  with  respect  to  the  value,  <h  ,  at  midspan  of  the  plate  so 

that  TrH  rH  1 

N     =    pv  ei(wt  +  ♦  '    \J0    PQU)cos  6  dz  +  i    JQ   roU)  sin  8  dz 

Assuming  application  of  the  lift  force  at  the  trailing  edge  the 
moment  on  the  plate  may  be  written  as 
M  =  2bN 

r     H  H  -J 

M     -     2pVbei(vt  +$  )      J     [,()(z)cos  9  dz  +  i       "po(z)sin9dz    I 


7. 

It  is  now  assumed  that 


p3? 

I  0   f 

where  P-.  =  circulation  in  the  vortex 

u*  =  the  longitudinal  fluctuating  velocity  component  as  measured  by- 
probe  B 

and  f   =  vortex  shedding  frequency. 

From  reference  (4) 

w 

Vo  — f. 

where  V  =  the  velocity  of  the  separation  streamline.  So  that  the  assumption 
here  is  that 

"I  ~vs 

Substituting  for  Pn(z)  the  moment  may  be  written 

H  MT 


rv  pVbe"~      '  '  \  J 

0 


i(wt  +4  )    [   "B      fl  A 

T   I     — =-  cos  C7  dz  +  1 


J0  ¥^^] 


Now  assuming  the  moment  coefficient  to  be  a  real  quantity  and  equating 
the  two  expressions  for  moment  we  have 

M 

_  .  I  •PKUU  I      \   1 


o 


[f bHvQ    \     u^5"   cos  9  dz  1 


and 


uj':  sin  9  dz  =  0 


Here  the  assumption  is  that  the  wake  is  symmetrical  behind  the  plate,  so 
that  the  integral  of  an  odd  function  is  zero.  Therefore,  the  forcing  moment 
coefficient  is  proportional  to  the  integral 


;. 


8. 


u*2  cos  8  dz 

0   * 

where  u*2"  is  the  longitudinal  fluctuating  velocity  component  measured  by 
the  traversing  turbulence  probe  (probe  B). 

B.  Turbulence  Probes 

Reference  (4)  describes  the  development  and  theory  of  the  total  head 
tube  in  detail.  Briefly  the  following  is  a  resume  of  the  theoretical  con- 
siderations involved  in  the  measurements. 

In  general  the  instantaneous  response  of  a  total  head  tube  may  be 

written 

H    =  p  ,  p'  ,  (u  +  u*)2  ,  (v  +  v')2  ,  (*  +  w')2 
tg     I        )(  2g         2g         2g 

where  p  is  the  pressure  and  u,  v,  and  w  are  the  velocity  components  in  the 
coordinate  direction.  The  barred  quantities  are  mean  values  and  the  primed 
symbols  are  fluctuating  values. 

From  the  geometry  of  the  probe,  it  is  a  good  assumption  that  the 
turbulence  probe  is  much  more  sensitive  to  the  velocity  components  parallel 
to  the  mean  flow  than  to  those  which  are  perpendicular  to  the  mean  flow 
(i.e.,  the  terms  vl2/2g  and  wt2/2g  are  neglected).  Also  low  frequency 
response  of  the  piezoelectric  crystal  makes  the  turbulence  probe  insensitive 
to  mean  quantities. 

Therefore  the  instantaneous  total  head  response  becomes  with  these 
assumptions: 

u    =  £*  +  2u"ut  +  HH 
tg     J     2g     2g 

The  D.C.  portion  of  this  signal  (the  time  average  of  the  fluctuation) 
is  given  by: 

Htg  "  "2g 


The  D.C.  level  is  removed  by  a  filter  so  that 

tg  "  ntg     |    g     2g  "  2g 

An  important  assumption  at  this  stage  is  that  all  terms  on  the  right 
hand  side  of  the  equation  can  be  ignored  with  respect  to  uuf/g»  whereupon 

tg      tg        g 

and  the  electrical  signal  as  measured  by  the  turbulence  probe  is  proportional 
to  the  product  of  uuf,  the  product  of  the  mean  longitudinal  wake  velocity  and 
the  longitudinal  fluctuating  velocity  component. 

C.  The  Wave  Analysis  System 

The  two  electrical  signals  measured  by  the  two  turbulence  probes  located 
in  the  wake  of  the  plate  and  positioned  along  the  span  were  analyzed  by  the 
wave  analysis  system. 

In  general  the  wave  analysis  system  accepts  one  or  two  periodic  or 
random  signals  as  inputs  and  provides  as  output  the  power  density  of  either 
signal  and  the  cross  power  density  expressed  as  a  complex  number.  The 
power  density  is  a  representation  of  the  amount  of  power  in  the  input  signal 
contained  in  a  given  bandwidth.  The  cross  power  density  is  a  measure  of  the 
power  (contained  in  a  frequency  bandwidth)  common  to  the  two  input  signals. 

If  the  probe  signals  were  harmonic  they  could  be  written  as  functions 

of  time  as  fA(t)  a^nci  fR(t)  and  their  transforms  as  F.Cnj)  and  FR(n2).  If 

F.(n1)  is  multiplied  by  FR(n2),  the  complex  conjugate  of  FR(n2),  and  the 

product  is  averaged  in  time,  only  products  of  sommon  frequencies  remain  so 

£  (n)   _-  --r     (n) 

that  ^  AB   =  F.(n)FB(n)  =  cross  power  density.  ^  ._   may  also  be 

written 


10. 


AB 
This  can  be  shown  since 


§  .„<n)    =    Vpsd.    VpsdI 


,4 


F,(n)     =    i      J         f,(t)  e-inW't  dt 


T-i 


and 


F2(n)     =    i       J"       f2(t)  e"1^^  dt 


^        j        f,(t)f2(t+y)dt     =    f f     J     Mt>*     ^    F2(n)ein0)i(t  +T) 
-Ti  -Ti 


4  -ij 


2* 


J    P,(d).**"i*     i        J        f,(t)elaa'lt  dt 


-T 


■•4 


therefore 


2* 

f,(t)fa(t+r)dt     =      ^^F7(n)F2(n)]einu>^ 


-Ti  >t-»  -«o 


■4 

and 

Ti 

I     J     .-^^drj    J      f,(t)fi(t  + 


F,(n)F2(n)     =    J  e_lna)'rdr|  f,(t)fa(t  +  T  )dt 


T 


-14  -x* 


but 
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1  rTi 


$    M   j    .-**<****   |    f   f,(t)f2(t  +  7)dt 


"T£  '     "Ti 


and  since  machine  operation  gives 
F,(n)  =  V  PSDA 


then 


F2(n)  =  n/pSDq 


4 


AB  A    V    B 

The  TP  625  wave  analyzer  system  accomplishes  this  by  passing  the  random 
input  signals  fA(t)  and  fn(t)  through  a  set  of  amplitude  and  phase  matched 

A  D 

band  pass  filters  before  multiplication.  In  order  to  preserve  the  phase 
angle  ©  between  the  two  input  signals,  two  modes  of  operation  of  the 
multiplier  yield  the  real  and  imaginary  parts  of  a  complex  number  called 
the  co- spectrum  and  the  quad-spectrum  respectively  so  that 

n  £ 

1  J   \  -   l  ( co- spectrum)2  +  (quad-spectrum)2\ 
and 


Q       -1  quad- spectrum 
co-spectrum 

For  random  signals  these  quantities  have  the  units  of  power  density  or  power 

divided  by  the  filter  band-pass  at  the  center  frequency  of  the  filter. 

Now  if  it  is  assumed  that  the  signal  from  probe  B  is  correlated 

(i.e.,  similar)  with  the  signal  from  probe  A,  the  fixed  probe,  we  can  write 

u«  =  K(z)uJ 

where  K(z)  represents  the  proportionality  between  the  signals  as  a  function 
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of  their  spanwise  separation.  As  noticed  previously  the  output  of  the 


probes  are    u*u  so  thatj^  J  *\j   u*u»u2.  Note  that  it  is  tacitly 
assumed  that  u.  =  u„.  Therefore 


I*  J 


*J      U.fU*U2  uJu*U2 

A  a  Ad 


and  if  we  divide  by  I  rp  1  u2  we  have 


IUJ 


\u 

rv 

u'u'u2 
A  B 

UAUBu2 

1  AA 

u'u'u2 
A  A 

'     u2 

Substituting  for  u' 

fy  KUTu]2"     ■     K^Tu]2" 


AB 


Iaa"2  S^-3? 


~    s? 


Since  wave  analyzer  operation  gives 


i 


=     PSDA 
AA  A 


then 

I    2 


II 


AB 


•N/ 


PSDA   u2  B 

A 

and  referring  to  the  discussion  of  the  forcing  moment  coefficient  where 
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H 


m.  /v 
o 


[_fbHv]  -1   1   ^  cos  9  dz 


•  \\.\ 


This  suggests  a  method  of  evaluating  the  integral.  |CP     and  PSD.  values 
are  available  as  outputs  of  the  wave  analysis  system  and  values  of  u  are 
easily  measured  in  the  tunnel.  Values  of  9  are  available  as  shown  from 

r\     .   -1  quad- spectrum 
co-spectrum 
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IV  PROCEDURE 

A.  Total  Head  Turbulence  Probe  Calibration 

The  total  head  turbulence  probes,  electrometer  preamplifiers,  D.C. 
filters,  and  the  turbulence  probe  signal  amplifiers  were  calibrated  as 
assembled  units.  The  two  system  outputs  were  adjusted  by  changing  the 
amplifier  gain  on  one  system  so  that  the  output  of  the  two  systems  were 
equal  for  identical  inputs.  Figures  5  and  6,  Appendix  A,  show  the 
amplified  turbulence  probe  signals  when  the  probes  were  subjected  to  a 
step  impulse.  See  reference  (2)  for  calibration  details. 

B.  Tunnel  Measurements 

The  fixed  turbulence  probe  was  initially  positioned  in  the  tunnel  at 
y/t  =  0  and  x/t  -  1.  See  Figure  15,  Appendix  A  for  description  of  coordinate 
system.  Alignment  of  the  traversing  turbulence  probe  (probe  B)  with  respect 
to  the  fixed  probe  (probe  A)  is  assured  by  the  method  of  attachment  of  the 
two  probes  and  the  proper  positioning  of  the  traversing  probe  seal  at  the 
bottom  of  the  test  section.  Measurements  were  taken  with  the  plate  free  to 
vibrate  at  ambient  velocities  of  8,  11,  and  14  fps,  and  clamped  in  a  fixed 
position  at  an  ambient  velocity  of  14  fps.  The  velocities  were  obtained 
from  the  contraction  meter  as  described  in  reference  (2).  Mean  wake 
velocities  were  measured  as  described  in  reference (3 )•  While  maintaining 
a  constant  ambient  velocity  by  operation  of  the  inlet  and  foot  valves,  the 
turbulence  probe  signals,  plate  vibration  signal  and  a  marker  signal  were 
recorded  on  four  channels  of  the  protable  tape  recorder  for  each  spanwise 
position  (z/t)  of  the  traversing  probe.  Prior  to  and  during  the  actual 
recording,  the  signals  were  monitored  on  oscilloscopes  to  insure  that  the 
proper  signals  were  recorded.  The  pressure  in  the  tunnel  was  maintained 
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at  the  same  value  by  operation  of  the  foot  valve  and  was  kept  sufficiently 
high  (3  psi)  to  prevent  cavitation. 

C.  Data  Analysis 

The  recorded  signals  were  analyzed  with  the  Technical  Products  Company 
TP  625  Wave  Analyzer  System.  Each  run  consisted  of  a  recording  of  the  four 
signals  (two  fluctuating  total  heads,  plate  vibrational  amplitude,  marker) 
at  a  particular  z/t  and  velocity.  For  each  run  the  probe  signals  were 
rerecorded  on  a  Sanborn-Ampex  2000  continuous  loop  tape  recorder  for  con- 
tinuous play  back  into  the  wave  analyzer  system.   Prior  to  each  analyzation, 
the  wave  analyzer  system  was  carefully  calibrated  with  a  100  cps,  O.46  volt 
peak  to  peak  oscillator  signal.  The  gains  of  the  two  wave  analyzers  were 
set  the  same  and  the  oscillator  (TP  626)  output  indicator  was  adjusted  with 
a  reference  oscillator  to  ensure  accurate  frequency  readings.  The  gains 
of  the  two  analyzers  TP  627  were  matched  and  the  Power  Spectral  Density  set- 
tings of  the  multiplier  TP  645  were  adjusted  to  give  the  same  output  for 
identical  inputs.   In  addition  the  phasing  between  the  identical  input 
signals  was  zeroed  using  the  phase  zero  dial.  Calibration  in  this  manner 
ensures  that  the  relative  power  levels  in  the  signals  are  preserved  in  the 
output,  and  that  zero  phase  angle  is  indicated  for  identical  input  signals. 
The  output  of  the  wave  analysis  system  is  in  volts  and  is  only  a  relative 
measure.  A  50  cps.  band  width  filter  and  a  long  averaging  time  (2  seconds) 
were  chosen  for  the  analyzation.  The  filter  and  2  second  averaging  time 
were  selected  to  give  smooth  curves  from  which  data  points  could  be 
accurately  determined. 

For  z/t  =  1  and  9  at  each  velocity,  the  power  spectrum  density  for 
each  signal  and  the  co-spectrum  and  quad-spectrum  were  obtained.  For  the 
remaining  valves  of  z/t,  the  power  spectrum  density  of  "Bw  (the  movable 
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turbulence  probe),  co-spectrum  and  quad- spectrum  were  obtained. 

The  horizontal  scale  of  the  x-y  plotter  was  adjusted  so  that  a  full 
sweep  covered  from  50  to  200  cps.  The  vertical  scale  of  the  x-y  plotter 
was  set  to  1  volt  per  inch.  Also  for  z/t  =  9  at  each  velocity,  the  power 
spectrum  density  "B"  was  obtained  using  a  1  cps  band  width  filter  and  an 
averaging  time  of  0.2  seconds  to  indicate  the  exact  frequency  of  the 
energy  in  the  wake.  Similarly  for  the  plate  motion  signal,  a  spectrum  was 
obtained  to  pinpoint  the  frequency  of  plate  vibration. 
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V  PRESENTATION  AND  DISCUSSION  OF  RESULTS 

The  outputs  of  the  Wave  Analyzer  system  are  the  power  density  spectrum 
of  a  single  signal,  the  co-spectrum,  and  the  quad-spectrum  of  the  two  signals 
taken  jointly.  These  are  presented  for  each  run  as  Appendix  C» 

For  analyzation  purposes,  frequencies  of  75  cps  for  runs  with  the  test 
plate  free  to  vibrate  and  10$  cps  for  runs  with  the  test  plate  fixed  were 
chosen.  These  frequencies  were  selected  because  they  were  the  respective 
points  of  peak  power  as  indicated  by  the  plots  of  power  spectrum  density 
(PSD)  and  hence  represent  the  frequency  of  vortex  shedding.  This  reason 
for  frequency  selection  is  reinforced  by  the  PSD  of  the  vibrating  plate 
signals  which  also  has  its  peak  power  at  75  cps.  See  the  power  spectrum 
density  for  vibrating  plate  signals,  runs  11,  16,  and  21,  Appendix  C. 

Power  density  spectra  were  obtained  from  the  wave  analyzer  system  with 
narrow  filter  settings  as  well  as  with  wide  filter  settings  for  runs  6,  11, 
16,  and  21.  These  narrow  filter  power  density  spectra  show  more  clearly 
the  exact  frequency  of  peak  power  than  do  the  power  density  spectra  made 
with  wide  filter  settings. 

Table  1,  Appendix  B,  lists  the  values  of  PSD,  co- spectra  and  quad- 
spectra  used  in  the  following  analysis. 

Power  density  spectra  for  probe  A  (fixed  probe)  are  shown  in  Appendix  C 
for  runs  1,  6,  7,  11,  12,  16,  17,  and  21.  These  are  all  runs  with  probe 
separation  (z/t)=  1  and  9«  It  can  be  seen  that  for  each  velocity  PSD.  is 
the  same  for  all  probe  separations  z/t.  This  indicates  that  at  small  probe 
separation,  one  probe  does  not  affect  the  signal  of  the  other  probe  and 
that  the  signal  from  the  turbulence  probe  A  is  stationary. 

The  power  density  for  the  fixed  probe  A  increases  in  magnitude  at  the 
plate  vibration  frequency  with  an  increase  in  ambient  velocity.  However 
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the  power  densities  for  the  fixed  plate  (runs  1  through  6  of  Appendix  C, 
ambient  velocity  -  14  ft/sec)  are  higher  than  are  the  power  densities  for 
the  vibrating  plate  runs  with  the  same  ambient  velocity  as  the  fixed  plate 
runs  (see  runs  17  through  21  of  Appendix  C  for  vibrating  plate  runs  with  an 
ambient  velocity  of  14  ft/sec).  Thus  power  density  is  not  a  function  of 
ambient  velocity  alone. 

If  one  computes  the  self  excitation  parameter  of  reference  (2), 
V#  =  — ^fro  t   for  each  run  made  in  this  study  and  compares  q>  ain^   V*» 
it  can  be  seen  that  <D   decreases  as  the  self  excitation  parameter  V* 
increases.  V*  is  the  ratio  of  the  transverse  velocity  of  the  trailing 
edge  of  the  test  plate  to  the  ambient  water  velocity. 

It  has  been  shown  in  Section  III  that  the  power  density  spectrum,  U>   , 
is  proportional  to  (u*u)2.  Thus  $.A/u2  is  proportional  to  (u*)2.  The 
quantity,  (uf)2/f,  may  be  considered  a  measure  of  the  strength  of  the 
individual  vorticies.  Thus  ?  7u2f  is  proportional  to  the  strength  of 
vortices  being  shed.  Figure  19  shows  the  strength  of  the  vortices 

shows  that  the  vortex  strength  drops  off  rapidly  with  an  increase  in  the 

self  excitation  parameter  V*.  Thus  the  increase  in  forcing  moment  co- . 

efficient  m.  with  an  increase  in  the  self  excitation  parameter  V*  cannot 

o 
be  explained  by  an  increase  in  the  strength  of  the  vortices  being  shed. 

For  the  case  of  the  fixed  plate  (V*  =  0)  it  appears  that  the  values 
of  vortex  strength  and  coherency  may  vary  with  the  ambient  velocity. 

The  power  density  spectra  for  the  signal  from  probe  B  at  small 
separations  (runs  1,  7,  12,  and  17  of  Appendix  C)  are  essentially  the  same 
as  the  power  density  spectra  for  the  signals  from  probe  A. 


, -2f  being  shed  plotted  vs.  the  self  excitation  parameter  V*.  Figure  19 
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As  shown  in  Section  III: 

\  Q)  .Rl  =  \  (co-spectrum)2  +  (quad-spectrum)2 1 

p   _  .   -1  quad-spectrum 
co-spectrum 

Table  1  lists  the  values  of  \<T>  ARI   and  9  f°r  all  runs.  As  explained  in 
Section  IV,  Procedure,  the  wave  analyzer  system  was  adjusted  so  that  the 
outputs  were  all  relative  to  the  100  cycle  sine  wave  calibration  signal. 

The  Technical  Products  Company  wave  analyzer  system  when  in  the  power 
spectrum  density  mode  of  operation,  multiplies  the  input  signal  by  itself 
and  in  effect  becomes  a  squarer.  Time  averaging  then  eliminates  all 
products  except  those  of  common  frequencies.  Thus  the  power  spectrum 
density  of  signal  A  =§AA'V  (uJu)2*  likewise  $  BB  ^  (uju)2.  When  the 
wave  analyzer  system  is  used  to  obtain  the  cross  power  spectrum  density 
the  two  input  signals  are  again  multiplied  and  time  averaged. 

When  interpreting  the  machine  output  we  must  remember  that  all  <J>  fs 
are  averages  over  a  finite  bandwidth.  Thus  when  we  look  at  <I>AA  or  5)BB 
with  the  filter  centered  at  frequency  n  ,  F  (n  )  or  FR(n  )  may  be  zero  yet 
Taa  or  trr  can  De  non  zero  due  t°  contributions  at  other  frequencies 
within  the  filter  pass-band.  When  looking  at  4>AB  *re  are  most  interested 
in  its  value  at  the  vortex  shedding  frequency.  If  this  frequency  is  only 
slightly  different  at  the  locations  of  probes  A  and  B,  and  the  random 
turbulent  background  at  surrounding  frequencies  is  negligibly  small,  (pAR 
can  be  zero  and  yet  H>AA  and  <PRR  would  be  non-zero! 

Since  J^AB|  /v/  (u^u2),  ^M  ~   (u^u)2  and  <^  RR  ~s  (u£u)2,  then 

\$  AB\~  f$  ,ul  *  \$BB\   ere  UA  and  UB  haVe  the  Same  frec*uencies* 
Therefore  \<£R|   is  a  measure  of  overall  signal  strength  as  well  as  a 

measure  of  the  similarity  of  the  two  signal  amplitudes.  To  eliminate  the 
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effect  of  overall  signal  strength  all  values  of  |q)  AR\  for  each  velocity 
can  be  normalized  by  their  respective  q)   values  (for  each  velocity  u>  . . 
is  stationary  as  shown  in  Appendix  C). 

Therefore:  ^>     1^ar|/^AA  an<^  "*"s  on^y  a  measure  of  the 
similarity  of  signal  A  and  signal  B.  Values  of  §)   are  tabulated  in 
Table  1,  Appendix  B.  Figure  17  is  a  plot  of  <£  vs.  z/t  for  different 
values  of  the  dimensionless  self  excitation  parameter  V*  =  - — • 
Figure  18  is  a  plot  of  Q   vs.  z/t  for  different  values  of  V*. 

Vibrating  plate  runs  at  an  ambient  velocity  of  1U   fps  were  repeated 
after  a  two  week  interval  to  investigate  the  reproducability  of  the 
experimental  data.  Figure  21  compares  values  of  a>  and  Q  vs.  z/t.  It 
can  be  seen  from  this  figure  that  the  cross  power  spectra  compare  very 
favorably,  but  a  phase  angle  variation  is  to  be  expected. 

From  Figure  17  it  can  be  seen  that  the  correlation  of  the  turbulence 
probe  signals  for  vibrating  plate  runs  decreases  as  probe  separation 
increases  which  is  what  one  would  expect.  Figure  17  also  shows  that  the 
normalized  cross  power  density  increases  for  vibrating  plate  runs  with  an 

increase  in  the  self  excitation  parameter  V*. 

x  * 

Figure  19  shows  that  the  normalized  cross  power  density  <T)   increases 

as  the  self  excitation  parameter  V*  increases. 

Figure  18  shows  that  for  z/t  <  5  for  vibrating  plate  runs  there  is  a 
gradual  increase  in  0  as  z/t  increases,  but  there  is  no  significant  change 
in  8  for  an  increase  in  the  self  excitation  parameter  V*.  For  probe 
separation  z/t  >  5  for  vibrating  plate  runs  9  increases  with  an  increase 
in  z/t  and  decreases  as  the  self  excitation  parameter  V*  increases. 

Figures  17,  18,  and  19  show  that  the  shedding  of  the  vortices  becomes 
more  correlated  and  more  in  phase  as  the  self  excitation  parameter  V* 
increases. 
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For  the  self  excitation  parameter  V*  =  0  (fixed  plate)  with  a  velocity 
of  LV  ft/sec,  Figure  17  shows  ©   to  decrease  to  near  zero  at  z/t  j^  7  and 
&  to  increase  for  values  of  z/t  >  7.  For  V*  =  0  Figure  18  shows  ©  to 
increase  gradually  for  values  of  z/t  4*   6  and  9  to  increase  rapidly  to 
values  >  180°  for  values  of  z/t  >  6.  This  indicates  that  for  V*  =  0 
there  may  be  a  spanwise  cellular  structure  with  alternate  cells  180°  out 
of  phase.  Assuming  symmetry  the  wavelength  of  the  cell  would  be  14.  body 
thicknesses.  Roshko  (6)  reports  a  18  diameter  wavelength  for  measurements 
behind  a  circular  cylinder  at  a  body  Reynolds  number  of  80. 

As  shown  in  Section  III: 

■V  ~  (bHrvv8)-1      'f  AB*   cos  9  dz 

Values  of  (Id)  A   cos  8  )/(  4)    )  were  plotted  vs.  z/t  and  inte- 
grated from  z/t  »  1  to  z/t  =  9.  It  was  assumed  that  the  value  of  this 
integral  is  proportional  to  the  integral  evaluated  over  the  entire  plate 
span. 

m_  was  calculated  by 

m„  ^  (fVu2)    /    -— >—   coscJd(z/t) 


f 
o 


1    £AA 


These  values  of  mf  for  the  vibrating  plate  runs  are  plotted  vs.  V*  in 

o 
Figure  20,  Appendix  3.  The  dotted  line  in  Figure  20  is  a  plot  of  m_  vs. 

o 
V*  as  determined  from  reference  (2).  It  can  be  seen  that  these  two  plots 

of  m_  are  in  qualitative  agreement.  In  an  effort  to  determine  any  pre- 

dominence  of  either  1<£  Anl2/  $aa  or  9  in  the  integrand  of  the  equation 

used  to  determine  rar  (i.e.,  is  the  amplitude  coherence  or  the  phase  eo- 

o 
herence,  the  predominant  factor  in 
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Jlrfi     \* 
■J  '     cos  Q  d(z/t)?  )  it  was  assumed  that  cos   9    was 

constant,  hence 


ZA=9     HJ 


mf    ^    (fVu*)'1  J 


The  resulting  nu  plotted  vs.  V*  is  shown  in  Figure  20. 


Next  it  was  assumed  that  (  \4?  \  )/(  $   )  was  constant,  hence 
m  r^j    (fVu2)   I  cos  9  d(z/t) 


The  resulting  m_  plotted  vs.  V*  is  shown  in  Figure  20. 
o 

The  conclusion  drawn  from  this  fugure  is  that  it  is  impossible  to 
determine  the  predominate  factor  in 

m,  „    (fVu*)"1    -iil—  cos8d(z/t) 
f°  J  f  AA 

Undoubtedly  both  (rrAR|  )/(4!A.)  and  cos  9  are  important.  Indeed  it  may 

be  said  that  amplitude  correlation  and  phase  correlation  have  no  significance 

unless  presented  together.  Phase  correlation  has  no  significance  unless 

there  is  some  amplitude  correlation  and  vice  versa. 

The  calculated  values  of  m   for  the  fixed  plate  runs  assuming  that 

o 
the  integrals  evaluated  from  z/t  =  1  to  z/t  =  9  were  proportional  to  the 

same  integrals  evaluated  over  the  entire  plate  were  not  considered  valid. 

However,  by  extrapolation  of  the  integrand  over  the  entire  length,  assuming 

a  cellular  structure  with  alternate  cells  180°  out  of  phase,  values  of  m. 

o 
obtained  are  in  qualitative  agreement  with  reference  (2)» 

The  co-spectrum  of  the  vibrating  plate  runs,  runs  7  -  21,  Appendix  C, 

show  that  some  correlation  exists  at  li*.0  cycles.  This  correlation  was  not 
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investigated,  but  it  could  be  the  correlation  of  the  first  harmonic  of  the 
vortex  shedding  frequency.  The  vibrating  plate  spectra  for  runs  |k  2.1 
Appendix  C,  show  power  at  about  1A0  cycles,  and  this  could  explain  the 
correlation  at  1A.0  cycles. 
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VI  SUMMARY  AND  CONCLUSIONS 

Spanwise  cross  correlations  of  the  longitudinal  turbulent  velocity 
component  in  the  early  wake  of  a  flat  plate  were  measured  in  water  under 
conditions  of  flow-induced  torsional  vibration  of  the  plate  in  order  to 
ascertain  the  nature  of  the  known  dependence  of  the  forcing  moment  upon 
the  ratio  of  vibratory  plate  velocity  to  free-stream  fluid  velocity. 

The  following  conclusions  are  drawn: 

1.  The  increase  in  the  forcing  moment  coefficient  m_  with  the 

o 
self  excitation  parameter  V*  is  due  to  an  improvement  in  correlation 

or  coherency  of  the  wake  along  the  span  rather  than  to  increasing 

strength  of  the  individual  vortices  being  shed. 

2.  The  forcing  moment  coefficient  can  be  expressed  in  the  form 

—  1  f*\$S 

J0       $AA 

3«  Both  the  absolute  magnitude  of  the  cross  power  density  q> 

and  its  phase  angle  9  are  important  in  explaining  the  relationship 

between  the  forcing  moment  coefficient  m_  and  the  self  excitation 

o 
parameter.  Spanwise  variations  in  phase  angle  are  shown  to  be 

reduced  with  increasing  V*  and  coherence  of  vortex  strength  at  the 

vibrational  frequency  is  increased. 

k»     Measurements  of  the  spanwise  wake  structure  under  conditions 
in  which  plate  vibrations  are  prevented  are  consistent  with  a  cellular 
structure  with  a  phase  shift  of  180°  between  alternate  cells. 

5.  In  order  to  minimize  self-induced  vibration,  engineers  should 
design  structures  in  such  a  way  as  to  eliminate  wake  coherency.  This 
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might  be  accomplished  by  a  random  spanwise  variation  of  trailing  edge 
geometry,  or  by  attaching  to  the  trailing  edge  short  independent 
sections  which  are  free  to  rotate  about  the  trailing  edge. 
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VII  RECOMMENDATIONS 

1.  In  the  co-spectrum  of  probe  B  some  indication  of  correlation  was 
observed  at  1A.0  cps.  It  is  recommended  that  this  be  investigated. 

2.  It  is  recommended  that  measurements  such  as  taken  in  this  study 
be  taken  so  as  to  include  observations  of  the  wake  along  the  entire  span. 
This  will  necessitate  a  new  turbulence  probe  supporting  mechanism. 

3.  It  is  recommended  that  a  5,  10,  or  20  cycle  bandwidth  filter  be 
obtained  for  the  Technical  Products  Wave  Analyzer  System.  These  more 
narrow  filters  will  permit  a  more  accurate  determination  of  the  frequencies 
of  peak  power  in  both  the  power  spectrums  and  the  cross  power  spectrums. 
The  available  one  cycle  filter  proved  to  be  too  narrow  for  the  random 
signal  involved. 
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TITLES  OF  PICTURES  ON  THE  FOLLOWING  PAGE 

Figure  1:   Water  tunnel  test  section  showing  electrometers, 
amplifiers,  DC  filters,  and  strain  gage  controls. 

Figure  2:   Left:  Sanborn  2000  tape  recorder. 

Center  bottom:  Precision  Instrument  Company  portable 
tape  recorder. 

Right:  Technical  Products  wave  analyzer  system. 
Figure  ^:   Probes  and  supporting  structure. 
Figure  l±:       Probes  and  supporting  structure. 
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FIG.    5 
PROBE  A  SIGNAL 
CALIBRATION 


FIG,   6 
PROBE  B  SIGNAL 
CALIBRATION 


FIG.   7 
PROBE  A  SIGNAL 
VEL.=  0  FT. /SEC. 


FIG.   8 
PROBE  B  SIGNAL 
VEL.=  0  FT. /SEC. 
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FIG.   9 
PROBE  A  SIGNAL 
VEL.=  11  FT. /SEC. 


FIG.   10 
PROBE  B  SIGNAL 
VEL.=  11  FT./SEC. 


FIG.   11 
PLATE  SIGNAL 
VEL.=  0  FT. /SEC. 


FIG.   12 
PLATE  SIGNAL 
VEL.=  8  FT./SEC. 


FIG.   13 

PLATE  SIGNAL 

VEL.=  11  FT./SEC. 


FIG.   U 

PLATE  SIGNAL 

VEL.*  14  FT./SEC, 
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TURBUIXNCI  PROBB  A 


AMBIENT  WATER  v 


*-   X 


The   origin  of  the  coordinate  system  is  located  at  the 
center  of  the  test  section  at  the  trailing  edge  of  the  test 
plate. 


Figure  15:  Turbulence  probe  arrangement,  cross  section 
of  test  plate,  and  coordinate  system 
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SYSTEM 


FIG.  16        FLOCK   DIAGRAM  Or    INSTRVMEiNTATION  FOR     ONE  OF  THE  TrfO 
TURBULENCE  PROSES 
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Figure  18:  Correlation  phase  angle  vs.  non-dimensional  spanvd.se 
distance 
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Figure  19:  Normalized  cross  power  density  and  vortex 
strength  vs.  self  excitation  parameter 
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Figure  20:     Forcing  moment  coefficient  vs.   self  excitation 
parameter 
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Figure  21:  Reproduced  results;  normalized  cross  power  density 
and  correlation  phase  angle  vs.  non-dimensional 
spanwise  distance 
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Table  1  Measured  Quantities 


z/t       PSDA       co-      quad-    —    1Q3     - 
Run  dimensionless  ,  _   »„  spectrum  spectrum  >   , .  "   N   ,..,./   \  /-,  /   \ 
distance    (volts)*  (volts)*  (volts)*  (radians)  ("/sec)  (ft/sec) 


1 

2 
3 

1 
3 
5 

3.25 
n 

1.55 
1.50 
1.15 

0.12 
0.20 
0.33 

plate 
fixed 

10.85 
tt 

tt 

14 

tt 

tt 

4 

7 

tt 

0.45 

-0.10 

tt 

n 

5 

8 

tt 

-0.55 

-0.17 

n 

tt 

6 

9 

tt 

-1.25 

-0.73 

tt 

n 

7 

1 

0.55 

0.31 

0.07 

3.93 

6.04 

8 

8 

3 

ti 

0.22 

0.06 

tt 

tt 

tt 

9 

5 

n 

0.17 

0.05 

n 

tt 

tt 

10 

7 

n 

0.10 

0.05 

tt 

tt 

tt 

11 

9 

n 

0.05 

0.05 

tt 

tt 

tt 

12 

1 

1.63 

0.90 

0.15 

4.96 

9.50 

11 

13 

3 

n 

0.58 

0.30 

tt 

tt 

tt 

14 

5 

tt 

0.39 

0.25 

tt 

n 

tt 

15 

7 

it 

0.21 

0.20 

tt 

tt 

n 

16 

9 

n 

0.05 

0.20 

tt 

tt 

tt 

17 

1 

2.0 

0.83 

0.00 

4.63 

8.50 

14 

18 

3 

n 

0.53 

0.21 

tt 

tt 

it 

19 

5 

fi 

0.18 

0.24 

tt 

n 

tt 

20 

7 

tt 

0.05 

0.17 

tt 

tt 

tt 

21 

9 

tt 

-0.10 

0.08 

tt 

tt 

tt 

*  Values  measured  at  105  cps  for  fixed  plate  and  75  cps  for  vibrating  plate. 


Table  2  Derived  Quantities 
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Run 

V* 

dimensionless 
velocity 

$  AB 
(volts) 

e 

(degrees) 

dimensionless 
power  spectrum 

$             3 

i^1 

w 

ul4/* 

1 

plate 

1.55 

4.4 

0.48 

4.07 

6.260 

0.0198 

2 
3 

fixed 
V*  =  0 

1.52 
1.19 

7.5 
16.0 

0.47 
0.37 

3-98 
3.13 

6.000 
3.560 

tt 

it 

4 

0.46 

167.5 

0.14 

1.19 

-0.539 

it 

5 

0.58 

198.8 

0.18 

1.53 

-0.832 

it 

6 

1.45 

210.3 

0.45 

3.81 

-4.770 

« 

7 

0.077 

0.32 

12.7 

0.58 

15.90 

4.950 

0.0150 

8 

ti 

0.23 

15.3 

0.42 

11.50 

2.570 

ri 

9 

n 

0.18 

16.4 

0.32 

8.80 

1.550 

it 

10 

it 

0.11 

26.0 

0.20 

5.50 

0.546 

it 

11 

it 

0.07 

45.0 

0.13 

3.57 

0.173 

tt 

12 

0.071 

0.91 

9.5 

0.56 

6.23 

5.560 

0.0181 

13 

ii 

0.65 

27.4 

0.40 

4.45 

2.570 

it 

14 

it 

0.46 

32.7 

0.28 

3.11 

1.220 

it 

15 

ii 

0.29 

43.6 

0.18 

2.00 

0.417 

ii 

16 

ti 

0.19 

73.2 

0.11 

1.44 

0.077 

it 

17 

0.052 

0.83 

3.6 

0.42 

5.83 

4.780 

0.0278 

18 

tt 

0.57 

21.7 

0.29 

4.02 

2.100 

it 

19 

« 

0.30 

53.0 

0.15 

2.09 

0.375 

it 

20 

n 

0.18 

73.6 

0.09 

1.25 

0.063 

it 

21 

it 

0.14 

146.3 

0.08 

1.11 

-0.113 

it 
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APPENDIX  C 


POWER  DENSITY  SPECTRA 
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